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Based on the recent free electron model for cohesion in 
narrow metallic constrictions by Stafford et al., we calculate 
the quantum fluctuations in the cohesive force versus elonga- 
tion. The fluctuations are dominated by states near the Fermi 
energy, thus explaining their apparently universal magnitude 
of order ef/\f- We present numerical results for the force 
fluctuations in a simple geometry and show that they are well 
described by the contributions of a few classical periodic or- 
bits in the Balian-Bloch trace formula for the density of states 
of transverse motion. 
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The mechanical properties of atomic size constrictions 
between two reservoirs of standard metals, for example 
Au, have be en the subject of intensive research in the past 
few years.u □ In particular, simultaneous measurements, 
of the cohesive force and the electrical conductance!-!™ 
have shown a striking correlation between the mechani- 
cal and transport properties. In a regime where the co- 
hesive force is linear in the elongation ('elastic stage'), 
the electrical conductance exhibits plateaus similar to 
those found in two dimensional quagtum point contacts 
in semiconductor heterostructuresMD By contrast, both 
the conductance and the elastic force rapidly decrease 
in the narrow regime between two conductance plateaus 
('yielding n stage'). These observations were originally 
explainedo in terms of abrupt atomic rearrangements 
which appear with increasing elongation-and indeed clas- 
sical molecular dynamics simulationa3~E3 seem to sup- 
port this point of view. On the other hand, the appar- 
ent similarity between the conductance plateaus found 
here and in semiconductor quantum point contacts to- 
gether with the strong correlation between conductance 
and force suggests that elastic and yielding stages in the 
cohesive force may appear even for a smooth constriction 
geometry as a result of changing the number of discrete 
transverse modes for the electrons.U 

This suggestion has recently-, been developed by 
Stafford, Baeriswyl, and Burki.ta It is based on view- 
ing the transverse eigenstates of the conduction electrons 
as delocalised chemical bonds which provide both con- 
duction and cohesion. Remarkably, a corresponding free 
electron model qualitatively accounts both for the aver- 



age cohesive force and for the abrupt steps in the force of 
order Ef/^f which are caused by successively cutting off 
the discrete electronic modes in the constriction. Now it 
is evident that any purely electronic model is an ideal- 
ization of the true experimental situation. For instance, 
possible atomic rearrangements which repeal themselves 
through a hysteretic behavior of the forceQ are neglected. 
A realistic electronic description moreover has to account 
for the strong deviations from a naive conductance quan- 
tization picture arising from the rather small value of the 
Fermi wavelength which makes the electronic motion sen- 
sitive to defects even on an atomic scale. tU Nevertheless, 
the simple free electron model of Stafford et al. not only 
provides an intuitive explanation of the observed behav- 
ior but also gives a correct order of magnitude estimate 
of the relevant forces, in the same way perhaps as the 
measured bulk modulus of many metals is roughly de- 
termined by the grousd state Fermi gas pressure of the 
conduction electrons Ji-a 

In our present work we discuss the fluctuations in the 
cohesive force which arise from the discreteness of the 
electronic eigenstates fps motion in the direction trans- 
verse to the elongation.Ej As shown by Stafford et al., the 
electronic cohesive force F of a constriction of length L 
is obtained from the T = free electron grand canonical 
potential 
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simply via F = —dfl/dL. Here p±(E,z) is the differ- 
ential density of states (DOS) for the transverse motion 
of the electrons at a given cross section (z = const.). 
Within a semiclassical description, which is valid as long 
as the Fermi wavelength Af is much smaller than the con- 
striction width R, the DOS may be split into an average 
p±(E,z) and a fluctuating contribution 5p±_(E,z). The 
latter vanishes if averaged over an energy range much 
larger than the spacing A(z) between successive trans- 
verse eigenstates. Assuming the deformation occurs at 
constant total volume V, the average cohesive force F 
which is associated with thepaverage DOS has only two 
contributions for Af <C i?Jl3'E3 



(2) 



Here the dominant contribution is associated with the 
(electronic contribution to the) surface tension, giving a 
cohesive force proportional to the change in surface area 
S with elongation L. In addition there is a universal 
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contribution 4 -P 1 to the average cohesive force which is 

F □ 

completely independent of the geometrylila It arises from 

the topological term in the Weyl expansion for the inte- 
grated DOSl-3 and leads to a weakening of the cohesion 
compared to the macroscopic surface tension contribu- 
tion. As is evident from equation the average co- 
hesive force does not depend on the detailed geometry 
of the nanosize constriction which will presumably vary 
significantly between different realizations of an experi- 
ment. Regarding the force fluctuations SF which arise 
from the fluctuating partr-af the DOS, it was found in 
the numerical calculations^ that the precise form of SF 
is specific to the shape of the cross section but does not 
depend on the detailed form of the constriction radius 
R(z) versus z. Remarkably, the rms amplitude of these 
fluctuations turned out to have a universal magnitude 
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independent of geometry. In order to understand these 
observations, we note first that 5F is determined by the 
fluctuations in ±h p i DP S of transverse motion. Follow- 
ing recent worltE^HLla we may therefore assume a simple 
cylindrical form R(z) = const, for the constriction with 
a radius R which scales like L^ 1 / 2 in order to fulfill the 
constraint of constant volume. Since the discrete eigenen- 
ergies e v for transverse motion all scale like R~ 2 , we have 
de v jdL — e u /L. It is then easy to show from ([!]) that 
in this simple geometry the fluctuations of the cohesive 
force at fixed total volume are given by 
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Here 5p±(E) is the fluctuating transverse DOS at the 
narrowest point of the constriction, which also deter- 
mines its conductance G. Now, in contrast to G which 
is a property of the states directly at the Fermi en- 
ergy, the cohesive force obviously depends on all the 
states with energies between zero and Ef- The numerical 
calculationsEjilj, however, which show that the force os- 
cillations are directly correlated with the conductance, 
indicate that 5F is dominated by the states near the 
Fermi energy. Indeed, this observation can be understood 
easily from equation (^), at least on a qualitative level. 
Since the fluctuating DOS Sp±(E) is a rapidly oscillating 
function which varies on a scale of the order of the mean 
level spacing A <C ef, the contributions to SF\y in (Q) 
from energies between zero and close to Ef cancel. It is 
only in a small range of several level spacings near the 
upper integration limit, where the factor (1 — E/ef) 1 ^ 2 
changes rapidly on a scale on which the DOS varies. 
Therefore only the contribution to SF from a few states 
below the Fermi energy survives in As a result, 

rms(SF) is expected to be of order £f/Xf independent 
of the constriction radius R, a property which has been 



verified, numerically up to values kpR = 200 by Stafford 
et al.E3Ej By contrast, if all the states from zero up to 
Ef were to contribute to the force fluctuations, SF would 
scale as the fluctuations in the total number of states 
N(ef) below the Fermi energy. Assuming Poisson statis- 
tics, appropriate for a classically integrable transverse 
motionH, one has rmsN(e F ) = N 1/2 (ef) ~ k F R- The 
force fluctuations would thus increase with the constric- 
tion radius, in contradiction with the results of Stafford 
et al@ 

For a quantitative confirmation of our arguments 
above, we have calculated numerically the force fluctu- 
ations which follow from ([|) for a cylindrical constric- 
tion with radius R{L) — Rq\J Lq/ L. The correspond- 
ing results for SF in units of the fundamental force 
£f/Xf ~ 1 nN are shown in Fig. [l] for a wire which is 
stretched from kpR — 30 down to kpR — 2. Clearly, the 
force fluctuations are independent of kpR, with an aver- 
age magnitude rms(SF) = 0.6 Ef /Xf which is of the sana© 
order as the topological contribution to the force 
The apparently complicated dependence of the force fluc- 
tuations on the elongation can in fact be simply under- 
stood as arising from only a few classical periodic orbits 
of electrons in the assumed circular cross section of the 
wire. Indeed, as is well known, the fluctuations in the 
DOS around its average value can generally be repre- 
sented in terms of a sum over classical periodic orbits 
in a semiclassical approximation. For integrable systems 
such a connection between the quantum mechanical DOS 
and classical mechanics was first found by Balian and 
Bloch.E^I Specifically, for a circle the periodic orbits are 
regular polygons. They may be characterized by their 
number of vertices v and their winding number w. Ob- 
viously we have v > 2w. If there is a common divisor 
n between v and w, the orbit is an n-fold repetition of 
a primitive periodic orbit (see Fig. || for some elemen- 
tary examples). Introducing an angle <p vw = Ttw/v, the 
length of a periodic orbit is L vw — 2vR sin cj) vw from sim- 
ple geometry. The oscillating contribution toJiic DOS of 
a circular billiard can then be represented ascil 
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Here A = is the unit of energy and f vw = 1 for 

v = 2w or f vw — 2 for v > 2w the number of differ- 
ent periodic orbits through an arbitrary point within the 
circle. Using the semiclassical approximation (^) in our 
expression (||) for the force oscillations, we find that the 
details of the exact numerical result for SF are essentially 
explained by including only the three simplest periodic 
orbits v = 2, 3, 4, w = 1 in the circle (see Fig. l]). Ex- 
tending the series to all the 16 orbits with v < 10 and 
w < 2, the agreement between the semiclassical and fully 
quantum mechanical calculation becomes essentially ex- 
act. The fact that only a few periodic orbits are required 
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to describe the force oscillations is a consequence of the 
integration over all energies in (^), which suppresses the 
DOS fluctuations on very short scales associated with 
longer periodic orbits. 

Finally, it is interesting to point out that the approx- 
imation of simply adding the Weyl and trace formula 
contribution to the DOS, which apparently works very 
well for our present problem, is not always-valid as has 
been shown very recently by Bhaduri et alEB 

In conclusion, we have studied the quantum fluctu- 
ations in the cohesive force of metallic nanowires which 
arise from the discreteness of the electronic motion in the 
transverse direction. It has been shown that only a few 
states below the Fermi energy contribute to these fluctua- 
tions, supporting the prediction of universal force fluctu- 
ations rms(6F) = const. -ef/\f first found by Stafford et 
al.tj Unfortunately even in our simple geometry we have 
not found an analytical derivation of this result. Indeed, 
in our semiclassical approach, if we substitute equation 
(0) into (|j) and then change the integration variable to 
E/ A it is not evident that the fluctuations are indepen- 
dent of kpR. Therefore the issue of force fluctuations and 
in particular their geometry dependence deserves further 
investigation. The detailed structure of the force fluc- 
tuations reflects the classical periodic orbits in a given 
cross section which we have assumed to be circular here. 
It would clearly be of interest to generalize our results 
to the case of chaotic motion in the transverse direction 
where the eigenvalues obey random matrix theory rather 
than Poisson statistics. Similarly to the situation of per- 
sistent currents in ballistic billiard structures, we expect 
that the force fluctuations will be smaller if the trans- 
verse motion is chaoticn3 Experimentally, the observa- 
tion of the quantum fluctuations in the cohesive force as 
well as the closely related charge fluctuations of order e 
predicted very recently by Kassubek et al.t3 would con- 
stitute a crucial test of the electronic model of cohesion in 
metallic nanowires as opposed to a classical mechanical 
model. 

We acknowledge David A. Wharam for carefully read- 
ing the manuscript. 
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FIG. 1. Comparison of the fluctuating part of the cohe- 
sive force SF with SF SC , calculated with three periodic orbits 
(v = 2,3,4, w = 1). 
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FIG. 2. Closed classical periodic orbits in a circlular bil- 
liard with reflecting walls. Between consecutive reflections 
the trajectory follows straight lines. The winding number w 
and the number of vertices v of the particular orbits are given 
as tupel (v,w) (after Balian and BlochEj). 
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